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1 Introduction 


The dual realizations that D-branes are an intrinsic component of string the¬ 
ory and that the five distinct vacua of string theory are actually different 
limits of eleven-dimensional M-theory have had dramatic implications for 
cosmology. Both developments widen the range of options open to model- 
builders, but they also increase the number of constraints that a successful 
“string” cosmological model must satisfy. In particular, branes have given 
rise to a huge number of new models, which typically postulate that our 
apparently three dimensional universe is actually a 3-brane embedded in a 
higher dimensional background @ . Even in “standard” cosmological models, 
however, branes are expected to be in thermal equilibrium with other matter 
in the early universe, and must therefore be incorporated into any complete 
cosmological scenario. 

Brane gas cosmology is devoted to understanding the cosmological role 
of the “gas” of branes that inhabited the early universe. A particular focus 
of this program has been whether brane gas cosmologies single out a special 
number of “large” dimensions. In 1989, Brandenberger and Vafa argued that 
cosmological models containing a gas of winding strings can produce, at most, 
three macroscopic dimensions |Q. This work has been tested and extended 
in a number of ways since then BB. BBS 0.0. Q . The key observation 
underlying this work is that extended, p dimensional objects have p -|- 1 
dimensional worldvolumes, and from a purely topological perspective these 
worldvolumes will only intersect in 2(p + 1) spacetime dimensions or less. 
Thus, strings can only interact in three spatial dimensions (or less), whereas 
2-branes can find each other in up to hve spatial dimensions. One obviously 
hopes that Brandenberger and Vafa’s conclusion that strings naturally lead 
to a three dimensional universe is not undermined when branes are added to 
the picture. Indeed Alexander, Brandenberger and Easson have argued 
that when branes are included, strings will still dominate the evolution of the 
universe at late times, so the basic Brandenberger-Vafa mechanism survives. 

This paper extends brane gas cosmology in two important ways. First, 
rather than work with ten dimensional string theory, we take eleven-dimensional 
M-theory as our starting point. This is a natural basis for any theory that 
hopes to understand the origin of the three large dimensions, since string the¬ 
ory is itself a compactihcation of M-theory, with the compact direction pro¬ 
viding the string dilaton |0. Consequently, string theory implicitly assumes 
that one direction is already on a different footing from the rest, and in the 
long run we would hope to explain this rather than inject it as a hypothesis. 
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Second, rather than using solely thermodynamic arguments, we study the 
dynamical evolution of an M-theoretic universe containing M2-branes and a 
gas of supergravity particles. We ignore the M5-branes that are present in 
the M-theory spectrum, since they can annihilate in the full 11 dimensional 
M-theory spacetime. We are particularly interested in cases where the dif¬ 
ferent directions (assumed to be toroidally compactihed) are anisotropically 
wrapped by M2-branes. We develop the equations of motion for the general 
case in which each of the 10 spatial dimensions is distinguishable from the 
others, and describe the conhguration of wrapped branes and anti-branes 
with a “wrapping matrix”. Although we assume a toroidal topology, numer¬ 
ical and analytical results 0 indicate that wrapping dynamics should be 


insensitive to topology, and hence we believe our conclusions extend to more 
realistic compactihcations. 

Strings are obtained when one dimension of the M2-branes gets wrapped 
on the 11th dimension, which is taken to be much smaller than the Planck 
scale. Similarly, D-branes in string theories can be obtained as various con- 
hgurations of M2- and M5-branes. Thus we no longer need to consider hve 
separate brane gas cosmology scenarios, one for each string theory: our goal 
is to see if we can use M-theoretic models to make generic predictions about 
the overall form of the universe. 

Motivated by the original Brandenberger-Vafa scenario, we focus on the 
late time asymptotics of universes with a subspace that has no wrapping 
modes, presenting both an analytic discussion and specihc numerical solu¬ 
tions. Here we are primarily concerned with what is possible, rather than 
what is probable, and focus on the most interesting initial states that are 
permitted by general topological arguments. In a subsequent paper we will 
discuss the thermodynamics of a brane dominated universe in more detail, 
in order to establish the relative likelihood of different initial states. Conse¬ 
quently, we pay most attention to the following scenario. We imagine that 
thermal fluctuations drive subspaces of various dimensions to momentarily 
expand to a larger than average size. If the fluctuation involves a hve dimen¬ 
sional subspace, interesting nontrivial dynamics can ensue. Fully wrapped 
2-brane/anti-2-brane pairs will generically annihilate in hve space dimen¬ 
sions, thereby removing their restoring force. This leaves partially wrapped 
2-branes that have one dimension wrapped along the expanding hve-space 
and one dimension wrapped along one of the other, smaller dimensions. From 
an ehective hve-dimensional perspective, such branes appear to be wrapped 
one-dimensional objects - strings - and naively appear to impede further 
expansion. Yet, we imagine that subsequent thermal huctuations within this 
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five-dimensional environment will, once again, drive various subspaces to ex¬ 
pand further, and if such a subspace should happen to have three (or fewer) 
dimensions, the effective wrapped strings will be able to annihilate, thus 
allowing the subspace to expand without further constraint. We therefore 
see the possibility of generating an interesting hierarchy with a large three 
dimensional subspace emerging from an intermediate hve dimensional sub¬ 
space, with all remaining dimensions small One of our goals is to study 
this scenario in detail to determine if the dynamics leverages the topological 
reasoning we have used and drives a hierarchy between the dimensions. 

In Section 2 we derive the equations of motion and discuss some of their 
general properties. In section 3 we present numerical solutions to the equa¬ 
tions of motion, and in section 4 we study the behavior at late times ana¬ 
lytically. In section 5 we consider string winding, to make contact with the 
original scenario of Brandenberger and Vafa. In section 6 we discuss the 
implications of our results, and outline possibilities for future work. 

2 Homogeneous brane gas dynamics 

Our models are governed by M-theory which is well-described by eleven¬ 
dimensional supergravity when the radii and curvature scales are larger than 
the eleven-dimensional Planck length. This assumption will be justihed in 
hindsight, as the cosmologically relevant solutions we obtain contain growing 
radii. We consider two types of matter. The hrst is the massless supergravity 
spectrum consisting of 128 bosonic and 128 fermionic degrees of freedom - 
we ignore massive modes since we expect that these will quickly decay. The 
second ingredient is wrapped M2-branes. Although M-theory also contains 
M5-branes, these will annihilate quickly in 10-1-1 dimensions and thus should 
not be signihcant to the late-time behavior. We assume that the M2-branes 
are at rest and ignore fluctuations on the brane volume, as these effects 
are subleading in our large radii assumption. Finally, to make the analysis 
tractable, we assume that the particle and brane gases are homogeneous. 

Using the metric ansatz 


d 




( 1 ) 


2 = 1 

the non-vanishing Christoffel symbols are 




( 2 ) 
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while the non-vanishing components of the Einstein tensor are 



where there is no implied summation on i in the second line. 

We hrst introduce a gas of massless supergravity particles, with energy 
density ps and pressure ps- For simplicity we take the gas to be homogeneous 
and isotropic, with a perfect fluid stress tensor 

T^^ = diag(-ps,Ps,---,Ps)- (5) 

The equation of state appropriate for d spatial dimensions hxes ps = ^Ps, 
while covariant conservation of the stress tensor requires 

const. 

~ Y{d+i)/d ’ 

where the volume of the d dimensional torus is simply 

V = l[27rR,. (7) 

i 

The second source of stress-energy in our model universe is a gas of 2- 
branes, wrapped on the various cycles of the torus. These are characterized 
by a matrix of wrapping numbers Wj, where we take to represent the 
number of branes wrapped on the [ij) cycle, while represents the number 
of antibranes. A single M2-brane is described by the Nambu-Goto action]^ 

S = -Ta j d^C\/-det (8) 

gap = (9) 

where the brane tension Ta = hes ■ This leads to the stress tensor 

= -Ta j d^^5^\X - X{0)V^^9 9'^f^d^X^dpX ‘'. (10) 

®See also a similar discussion of the action of a p-brane in fl^. 
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For simplicity we assume that the branes are at rest, and ignore any possible 
excitations on the brane worldvolumes. Then, for a single brane wrapped 
on the (12) cycle and uniformly smeared over the transverse the stress 

tensor is 

With a matrix of wrapping numbers, the non-zero components of the brane 
gas stress tensor are 

T\ = (12) 

k^l 

T\ = -^(27r)2^i?fci?,(W* + iV*fc) (13) 

k^i 

with no sum on i in the second line. 

We now insert these two sources of energy-momentum into the right hand 
side of the d-dimensional Einstein eqnations, 

(14) 

where the gravitational coupling is 


IGttG = . 


( 15 ) 


The time-time eqnation can be solved for the energy density and hence pres¬ 
sure of the super gravity gas, 

= (16) 

^ k^l '■ k^l 

After a little re-arrangement, the space-space eqnations yield the following 
set of second-order differential eqnations for the radii. 


Ri 

Ri 


87rGT2 

V 


2d+l 
d{d — 1) 


(27r)2^i?fci?Wfc/- 

k^l 


(27r)^ ^ ^ RkRi {Nki Nik) 

k^i 



E 


^k^i 

^k 


( 17 ) 
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2.1 Some exact solutions 

With no wrapped branes (Wj = 0) the model reduces to an 11 dimensional 
radiation dominated universe, and it is easy to hnd a number of exact solu¬ 
tions. In particular, with the ansatz 


Ri{t) = Ci {t - (18) 

we hnd the following solutions to (p^Tl): 

• Flat spacetime, with ai = ■ ■ ■ = = 0, in which case ps will vanish. 

• Kasner solutions, with = 1- The energy density ps 

vanishes, so these are vacuum solutions. Note that the exponents a* 
must lie in the interval [—1,1], and that at least one of the ttj’s must 
be negative. 

• Radiation dominated, with ai = ■■■ = ad = The energy density 
of the radiation is non-zero, hxed by the constraint (0). 

One can also hnd a solution with brane wrapping but no gravitons, in which 
we set Nij = n for all i j but take ps = 0. Using the same ansatz (|T8|), but 
setting all the Ri equal (so that q = c and «* = a), demanding that ps = 0 
implies 

op- _ 87iGT2n 

2(t-toy “ (Mp^‘ ^ ^ 

This equation is solved for all times by setting a = and choosing the 
value of c appropriately. Given these values, one can check that the Ri 
equations (|I^ are also satished. 

The ansatz ([I^) doesn’t give the most general solution to the equations 
of motion, since we ought to be able to specify 2d independent initial radii 
and velocities. 

2.2 General features 

To study the solutions of the held equations o, it is convenient to work 
in terms of new variables Xpt) = log (27ri?j(f)) introduced in [^. The held 
equations become 

A) + = SttG Qps + ^^(W, + . (20) 
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The volume of the spatial torus is 

V = , 

the energy density from brane tension is 




( 21 ) 


( 22 ) 




and the energy density from supergravity excitations is hxed by the constraint 

(23) 


0 ) 


1 


SirGps = -iih “ SitGpb- 


i¥=j 


As noted in P, this system of equations can be regarded as describing a 
non-relativistic particle moving in d dimensions. The particle has a coefficient 
of friction, given by V/V, due to the expansion of the universe. A position- 
dependent force acts on the particle, given by the right hand side of (^Ol). 
This force consists of two terms. The first term. 


Fy=8.G(-ps + ,_i 


Pb 


(24) 


is positive definite and is the same for every value of i. It drives a uniform 
expansion of the universe. The second term. 


= 


SttGT 


V 




(25) 


is either zero or negative and it suppresses the growth of dimensions with 
large wrapping numbers. In the M-theory context, this is the mechanism by 
which anisotropic wrapping numbers lead to an anisotropic expansion of the 
universe. Another way of stating this is to note that the equations of motion 
imply that 


V . 

Ai - Xj + —(Ai - Xj) = 8nG{pi - pj) 


(26) 


where pi is the pressure exerted on the dimension. Thus differential 
pressures, of the sort exerted by an anisotropic brane gas, lead to differential 
expansion rates. 
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Let us comment on some general properties of these equations.^ First 
note that the volume of the torus increases monotonically with time. This 
follows from rewriting the constraint (^) in the form 

+ IGttG (ps + Pb) ■ (27) 



The right hand side is positive dehnite for all nontrivial models, vanishing 
only for the trivial case of periodically identihed Minkowski space. We can 
choose the direction of time to make > 0. 

The rest of this paper is primarily concerned with models where the 
wrapping matrix Nij is anisotropic. Motivated by our discussion in section 
1, we assume that the spatial dimensions fall into three classes. We refer 
to a direction i such that Nij = Nji = 0 for all j as unwrapped. Directions 
i for which Nij and Nji are nonzero except for those j corresponding to an 
unwrapped direction are referred to as fully wrapped. Directions i where some 
of the Nij or Nji are zero for values of j which are not unwrapped are referred 
to as partially wrapped. That is, listing dimensions in the order unwrapped, 
partially wrapped, fully wrapped, the wrapping matricies we will consider 
look like 


/ ■ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

iV4,6 

W,7 

-W,8 


-W,10 

0 

0 

0 

0 


iV5,6 

iV5,7 

iV5,8 

N5,9 

iV5,10 

0 

0 

0 

iV6,4 

iV6,5 



-W.S 

7^6,9 

-W.IO 

0 

0 

0 

iV7,4 

iV7,5 

7^7,6 


-A^7,8 

^7,9 

-A^7,10 

0 

0 

0 

Nsa 

iV8,5 


Ngj 


Ns,9 

-As,10 

0 

0 

0 

iV9,4 

iV9,5 


Ndj 

^9,8 


-Ag.lO 

VO 

0 

0 

^10,4 

Nlop 

Niofi 

-A^10,7 

-A^10,8 

Nw,9 



The diagonal entries in the matrix are irrelevant. In this example the direc¬ 
tions 1 through 3 are unwrapped, 4 and 5 partially wrapped, and 6 through 
10 fully wrapped (assuming all the Nij written out explicitly in the above 
matrix are non-zero). 


®These results overlap with the work of Don Marolf |Q. 





Now consider what happens when m dimensions are nnwrapped and d—m 
dimensions are either partially or fnlly wrapped. We introdnce 

m 

\i = log (volnme of nnwrapped torns), (29) 

i=l 
d 

\i = log (volnme of wrapped torus). (30) 

Summing over the appropriate values of i, and using the dehnition of pb, we 
hnd the following differential equations for p and A: 

ii+(fi + A)fi = *’^'3 j (31) 

A + (A + A)A = (32) 

The overall character of the dynamics depends on m. For small m, both terms 
in the right hand side for A are positive, and if A is zero, then the second 
derivative must be positive, leading to a local minimum. Conversely, with a 
larger value of m, the right hand side of this equation has both a positive and 
a negative term, and both local maxima and minima are allowed. For p there 
can never be a negative term on the right hand side, so if these directions 
are initially expanding they will then expand forever. 

When d — 3m + 2 vanishes the brane tension does not contribute to the 
growth of the internal dimensions. This picks out a special dimensionality 
of the unwrapped subspace, for which the internal dimensions have nearly 
constant radii. When d = 10 this occurs for m = 4. 

Perhaps the most important feature of these equations, which we analyze 
in more detail in section is that an attractor mechanism governs the be¬ 
havior at late times. That is, the radii at late times are determined solely by 
the wrapping matrix, and are (up to some scaling symmetries discussed in 
section ^) independent of the initial radii and velocities. This is intuitively 
clear from examining the second term in the “force”, (^Sf) . Suppose one of 
the radii, say Ri, starts out with an unusually small (or large) value. This 
does not affect the force on Ri itself, since Ri appears in both the numer- 
ator and denominator of . But the force on the other radii will become 
larger (or smaller), in such a way as to restore a balance between the different 
radii. This suggests that the behavior at late times is solely determined by 
the wrapping matrix. 


A = 

A = 
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3 Numerical results 


We now present numerical solutions to the equations of motion (P!7[). The 
technique is to start with a set of initial radii and velocities, then evolve the 
system both forwards and backwards in time using a Runge-Kutta routine. 
Motivated by our discussion in section |I|, we assume a wrapping matrix 
with mi unwrapped dimensions, m 2 partially wrapped directions, and m 3 = 
d — mi — m 2 fully wrapped dimensions. For the specihc case mi = 3 and 
m 2 = 2 in a ten dimensional universe we have the wrapping matrix given by 
([28|). The wrapping numbers should be chosen based on an understanding 
of thermal fluctuations in the early universe. Pending this, for our numerical 
calculations we simply take the non-zero wrapping numbers to be randomly 
chosen integers Nij G {1,2,3}. We take the initial radii Ri to be randomly 
chosen between 1 and 5 Planck lengths, and we take the initial velocities Ri 
to be randomly chosen between -0.5 and -1-0.5. A typical numerical solution 
for 3 unwrapped dimensions, 2 partially wrapped dimensions and 5 fully 
wrapped dimensions, is shown in Fig. 1. A solution for mi = 3, m 2 = 4 and 
m 3 = 3 is shown in Fig. 2. 

The equations predict that the early universe is radiation-dominated, with 
the volume of the spatial torus being V ~ Rd./[d+i)^ Thus the universe orig¬ 
inates in a big bang, which we have taken to be at t = 0. The big bang 
is apparently highly anisotropic, with some radii increasing while others go 
to zero, but we do not trust our equations of motion at early times: the 
temperature diverges while many of the radii are sub-Planckian, so we do 
not expect supergravity to be valid. It would be interesting to understand 
what, if anything, resolves the initial singularity in M-theory. Presumably 
U-duality plays a role, perhaps along the lines that have been argued in string 
theory g. 

We do expect supergravity to hold at late times, due to the growing 
radii and falling temperature. The growing radii allow us to neglect brane- 
antibrane annihilation, which turns off as the universe expands, and excita¬ 
tions on the branes will get red-shifted away faster than the brane tension. 
We therefore trust our solutions to provide an accurate description of the sys¬ 
tem at late times, presumably until low-energy dynamics eventually stabilize 
the moduli. 

Thus the most signihcant feature of these solutions is that, after a rather 
long transition period, the unwrapped radii begin to grow as a universal 
power of t. The partially and fully wrapped radii also grow as universal pow¬ 
ers of t, but with different (smaller) exponents. This leads to the dynamical 
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Figure 1: Radii vs. time for three unwrapped dimensions (solid red curves), 
two partially wrapped dimensions (dotted green curves) and five fully 
wrapped dimensions (dashed blue curves). Distances and and times mea¬ 
sured in Planck units. 
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Figure 2: A plot with three unwrapped dimensions (solid red curves), four 
partially wrapped dimensions (dotted green curves) and three fully wrapped 
dimensions (dashed blue curves). 
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formation of a hierarchy between the wrapped and unwrapped directions. 
We will determine the power laws analytically in the next section. 

To conclude, let us comment on a less generic feature of these solutions. 
As we will see in the next section, the exponents appearing in the late-time 
power laws depend only on the values of mi, m 2 , m 3 . However the coeffi¬ 
cients in front of the power laws depend have complicated dependence on the 
wrapping numbers. A wrapping matrix of the form (1^) tends to make the 
partially wrapped dimensions bigger than the fully wrapped dimensions at 
late times. This effect can be seen in the hgures, and suggests that brane gas 
cosmology also leads to a sub-hierarchy in the sizes of the internal dimen¬ 
sions.!^ In the case of most interest, this suggests that one may end up with 
hve small dimensions, two intermediate-sized dimensions, and three decom- 
pactihed dimensions. Such a hierarchy has also been argued to emerge from 
weakly-coupled string theory [§. 

The intriguing possibility of the dynamical hierarchy comes with a num¬ 
ber of qualihcations, however. We need a better understanding of thermal 
fluctuations at early times, to determine what sort of wrapping numbers are 
reasonable. In addition, in this scenario some low-energy mechanism must 
eventually stabilize the moduli, and it is quite possible that this mechanism 
will erase any hierarchy that forms. 


4 Late time behavior 

To complement our numerical results, we proceed to study the behavior of 
the radii at late times analytically. Rather remarkably, there is an attractor 
mechanism at work, giving late-time behavior that is independent of the 
initial radii and velocities. 

We study the late-time behavior with mi unwrapped dimensions, m 2 
partially wrapped dimensions and m 3 fully wrapped dimensions. That is, we 
consider wrapping matrices of the form of (^). 

The late-time behavior of solutions to the equations of motion (pTj) is 

^Of course the unwrapped dimensions grow with a faster power law and always become 
much larger than the internal dimensions. 
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captured by the ansatz 

{ Qit^ i = 1,... ,mi (unwrapped) 

bit^ i = mi + 1,mi + m 2 (partially wrapped) 

Cit'^ i = mi + m 2 + 1,..., d (fully wrapped). 

(33) 

With this ansatz, the energy densities due to supergravity particles and brane 
tension are given by 


Ps 

Pb 


f const. \ 

T 2 


V V 


const. 


n “i n n ^^^miQ:+m2/3+m37 


UaiUhU (^^^mia+m2/3+m37 
X 




{d+l)/d 

(34) 


(35) 


There are several distinct classes of behavior at late times since, depending 
on the values of mi, m 2 , m 3 , certain contributions to the energy density 
become negligible at late times. The following four subsections exhaust the 
possibilities. 


4.1 Neglect ps 


If the universe expands rapidly enough then it is brane-dominated at late 
times, with 3> ps- We begin by analyzing this case, which turns out to 
be the one of greatest physical interest. 

Let us assume that every term in pe makes a comparable contribution to 
the energy density. This amounts to assuming that /3 = 7 . Substituting the 
ansatz (PBD into the equations of motion (P 0 [), and neglecting ps, the time 
dependence cancels out of the equations of motion as long as 

miCK + (m 2 + m 3 — 2)/3 = 2 . (36) 


The equation of motion for each unwrapped dimension becomes 

3 

a(miQ! + (m 2 + ms)P — 1) = - - MGt^p-Q . 

d — 1 


(37) 


On the other hand, summing the equations of motion over all partially and 
fully wrapped dimensions (that is, from i = mi + 1 to d) gives 


/d(m2 + m 3 )(mia + (m2 + m3)/d — 1) 


/ 3(m2 + m3) 

V d- 1 



STrGf^PB • (38) 
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Comparing ( 0 ) and (|HB|), one must have 


a 


{m2 + m^)(3 d — 3mi + 2 


The solution to (|36D and (^) is 


a = 


6{d — mi] 


13 = 1 = 


d{d — mi) + 4mi — 4 
2d — 6 mi + 4 


(39) 


(40) 


d{d — mi) + 4mi — 4 

Note that the exponents depend on the spatial dimension d and the number 
of unwrapped dimensions rui, but not on m 2 or m 3 separately. 

The late-time power law exponents are therefore hxed by d and mi. We 
now show that the coefficients in front of the power laws are also depend 
on the wrapping matrix. The ansatz (|55|) reduces the equations of motion 
to a set of algebraic equations. The equations of motion for the unwrapped 
dimensions are all identical, and give a single nontrivial equation. The equa¬ 
tions of motion for the partially and fully wrapped dimensions give another 
'm 2 + m 3 equations. In addition equation (^) must be imposed. Thus we 
have a total of m 2 -l- m 3 -|- 2 equations. Now let’s count the unknowns. Since 
tti appears only in the combination a*, there are a total of m 2 -|- m 3 -|- 3 
unknowns a, /? = 7 , H®*) Cj- Thus for generic wrapping numbers the 
equations are sufficient to determine all of the unknowns in terms of the 
wrapping matrix, up to a freedom to make rescalings 

which leave the equations invariant (for example, note that pe is invariant 
under this rescaling). 

This shows that there is an attractor mechanism at work. Modulo the 


scale symmetry (|4l|) , the late time behavior is determined solely by the wrap¬ 
ping matrix. That is, the system forgets about the initial values of the radii 
Ri and velocities Ri. 

Of course, this analysis assumes that we can ignore ps. This is only 
consistent if 

(mio; -|- m2/3 -f m^'f) {d + l)/d > 2 (42) 

or equivalently 

mi < 3d(d + 2)/(7d + 2). (43) 

This analysis also assumes that all terms in pe are comparable. By comparing 
to section |4.2|, this requires m 3 > m 2 . 
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4.2 Neglect ps and full-full wrapping 

For a highly anisotropic wrapping, the partially and fully wrapped dimensions 
will behave differently at late times. To analyze this possibility, assume that 
7 < / 5 . This makes the full-full wrapping terms in the brane energy density 
negligible at late times. Let’s also assume that ps is negligible. 

We get one equation by requiring that pe ~ namely 

mio; + (m 2 - l )/3 + (m 3 - 1)7 = 2 . (44) 

The equation of motion for the unwrapped dimensions reads 

3 

Q;(miQ; + m2/9 + msq — 1) = -- SnGt'^pB (45) 

a — 1 


while the sum of the equations of motion for the partially wrapped dimensions 
gives 


m2/3(miQ; + m2/5 + msy - 1) 


3m2 
d — 1 


SnGt'^PB 


(46) 


and for the fully wrapped 


m37(miQ; + m2/5 + m^'j - 1) 


3m3 
d — 1 


SirGt'^PB . 


( 47 ) 


This is sufficient to £x the exponents 

6m2m3 

G: = - 

(d - 4)17121713 + (d - l)(m2 + m3) 

^ _ 2 m 3 ( 3 m 2 -d+l) 

(d - 4)17121713 + (d - l)(m2 + m3) 

21712(31713 — d + 1) 

(d - 4)m2m3 + (d - l)(m2 + m3) 

A little equation counting shows that the coefficients of the power laws are 
determined by the wrapping matrix, up to the freedom to make rescalings 
by 

tti Qtti bi ^ 

and by 
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This analysis assumes that 7 < /3, or equivalently that m 2 > m 3 . It also 
assumes that ps is negligible, which is consistent if 

(mitt + 7712/3 + 77137) (c? + l)/(i > 2 ( 50 ) 

or equivalently 

m2m3 ^ d{d- 1) 
m2 + m3 7d + 2 

4.3 Neglect full-full wrapping 

Next assume that pg makes an important contribution at late times, but that 
7 < /3 so that the full-full wrapping terms in the brane energy density are 
negligible. 

Requiring that ps and pe scale like 1/t^ gives two conditions on the ex¬ 
ponents. Together with the equations of motion, this is sufficient to £x 

2 d{m2 + m3) + 277727773 

Q? =- 

d + 1 d{m2 -|- m3) — 477727773 

^ _ 2 7/(7772 - 27773) + 277727773 

d +1 7/(7772 -|- 7773) — 477727773 

2 7/(7773 — 27772) + 277727773 

d+1 7/(7772 -|- 7773) — 477727773 


The coefficients in front of the power laws are hxed up to the rescaling free¬ 
dom 




Q 


-mil 




(52) 


which leaves ps and the relevant terms in pB invariant. 

Neglecting the full-full wrapping is consistent if 7 < /3, or equivalently 
if 7772 > 7773 . By comparison with section [4.2| , we need 77727773/(7772 + m 3 ) < 
d{d — l)/{7d + 2) in order for pg to make an important contribution. 


4.4 All terms important 

Finally, consider the possibility that all terms in the energy density are im¬ 
portant. Plugging the ansatz (^) into the equations of motion the time de¬ 
pendence cancels provided that both pg and pe scale like 1/t^. If we assume 
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mi 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

a 



4 

7 

1 

5 

4 

23 

1 

7 

2 


5 

7 

13 

2 

11 

11 

77 

4 

33 

11 


1 

4 

1 

5 

1 

7 

1 

13 

0 

I 

II 

I 

II 

1 

11 

I 

II 

I 

II 

- 


Table 1: Exponents for d = 10 spatial dimensions as a function of the number 
of unwrapped dimensions mi, for universes with more fully wrapped than 
partially wrapped dimensions (m 3 > m 2 ). 


that all contributions to pe are equally important this hxes the exponents 


a = 


/3 = 7 = 


3d — mi 

mi(d + 1) 

1 


(63) 


d + 1 


Note that, just as in section |4.1|, the exponents depend on d and mi, but not 
on m 2 or m 3 separately. 

The equations of motion become a system of m 2 + m 3 + 1 algebraic equa¬ 
tions for the m 2 -|-m 3 -|-1 unknowns ])([ a*, c*. In this case there is a true at¬ 

tractor mechanism at work: the unknowns are uniquely determined since the 
scale symmetry ( 0 ) is broken. This analysis assumes that all contributions to 
the energy density are equally important at late times. By comparing to the 
other cases we analyzed, this requires m 3 > m 2 and mi > 3d(d-|-2)/(7d-|-2). 


To summarize, the late time behavior falls into four classes, depending on 
the values of mi, m 2 , m 3 . The various classes are characterized by whether ps, 
certain terms in pe, or both can be neglected at late times. The exponents for 
each of these cases are worked out above; the results for d = 10 and m 3 > m 2 
are given in Table 1. Moreover, we have shown that the coefficients in front of 
the power laws are (modulo possible scale symmetries) uniquely determined 
by the wrapping matrix. This means, for example, that the ratios of the sizes 
of the fully wrapped dimensions stay constant as the universe evolves. 

One important issue remains to be addressed - whether we expect our 
equations of motion to provide an accurate description of the dynamics at 
late times. Classical supergravity is only valid if all radii are larger than 
the 11-dimensional Planck length, so we should only trust solutions in which 
the radii are either constant or grow with time. For such solutions we are 
also justihed in neglecting massless excitations on the branes, which redshift 
away faster than the brane tension. We can also neglect brane-antibrane 
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annihilation (which will turn off as the transverse dimensions expand). The 
condition for growing radii is / 3,7 > 0. This restricts us to considering one 
of two cases: either the universes discussed in section |0 but with 


mi < 


d + 2 


or the universes discussed in section but with 

d — 1 

"^3 > • 


(54) 


(55) 


5 String gas cosmology 


Aside from the cases and (|55| ) some of the wrapped dimensions will 
shrink with time, and the supergravity approximation to M-theory eventually 
breaks down. To follow the evolution of the universe, one must take the U- 
duality group of M-theory on into account |K]. As an example of this 


analysis, we study the following rather special set-up, which will allow us 
to make contact with the original string gas scenario of ill- Consider a 
wrapping matrix in which the only non-zero matrix elements are iVj lo and 
A'lo,* for i = 1,..., 9. Regarding as the dilaton direction, this corresponds 
to type IIA string theory on T®, with a gas of fundamental strings wound 
around all directions of the torus. One can solve the equations of motion 
([TtP at late times with the following power-law ansatz for the radii. 


Ri{t) ^ i = 1,..., 9 

Rio{t) ~ R (56) 


This analysis was performed in detail in section [4.3| . One finds that ps and pe 
make comparable contributions to the energy density at late times. Moreover, 
one finds that Rio shrinks with time. This invalidates the use of an eleven¬ 
dimensional equation of state for the supergravity gas, which is implicitly 
encoded in the constraint (0), and which would have 


Ps ~ 



(57) 


Instead the gas is effectively ten dimensional, and we should take the energy 
density to scale as 


1 


( 58 ) 
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The energy from brane tension scales as 


~ ^{d-2)p • 

The equations of motion (pOf) require that both ps and pe scale like 1/t^, 
which hxes the exponents^ 

(3 = 1/A, 7 = - 1 / 2 . ( 60 ) 


This determines the late-time behavior of the M-theory metric. But since 
/?io is shrinking, we should really reinterpret this as a llA string solution, 
using 0 

C^4-theory = e"^'^/^d4rmg + (61) 

where 0 is the dilaton and ds^mg is the string-frame metric. Thus 


~ r^/^ (62) 

during ~ const. (63) 

This makes contact with the scenario of ii, where a gas of string winding 
modes leads to dimensions whose size (as measured in the string frame) is 
stabilized at around the string scale. Note that this scenario, which proposed 
string winding as a way to stabilize radii, in fact has a growing volume with 
respect to the M-theory frame. 


6 Conclusions 

We have presented solutions to the equations of motion of M-theory in which 
a hierarchy arises between the wrapped and unwrapped dimensions. More¬ 
over, these solutions can potentially differentiate between the wrapped di¬ 
mensions themselves, depending on the specihc form of the wrapping matrix, 
producing a “sub-hierarchy” of compact dimensions. While these results are 
very suggestive, we have left a number of open issues. 

First, while the wrapped directions grow far more slowly than the un¬ 
wrapped directions, they are not actually stabilized in the most interesting 
case with three unwrapped (and therefore macroscopic) directions. However, 

® Although it may seem we are mixing 11-dimensional equations of motion with a 10- 
dimensional equation of state, we are secretly using the fact that IIA supergravity can be 
obtained by dimensional reduction from M-theory. 
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the fluctuations envisioned by the Brandenberger-Vafa mechanism in a string 
dominated nniverse are actnally more likely to generate one or two macro¬ 
scopic directions, since these flnctnations are smaller and thus favored over 
those that remove the winding strings from a three dimensional snbspace. 
Brandenberger and Vafa resolve this problem by pointing ont that anthropic 
arguments make it very nnlikely that a nniverse with only one or two large 
dimensions conld actnally be “observed”. In the case discussed here, though, 
taking three completely nnwrapped dimensions results in a set of compact 
directions that expands more slowly than the one or two dimensional cases. 
Conseqnently, it is presnmably easier for the three dimensional models to be 
stabilized by non-pertnrbative effects not acconnted for by the snpergravity 
action we are working with. Thus it is at least possible that a three dimen¬ 
sional nniverse could eventually be picked out as a special case where one 
had both a permitted nnmber of “large” dimensions, and a stabilized set 
of compact directions, which is vital if we are to avoid constraints on the 
time dependence of Newton’s constant in the “effective” three dimensional 
nniverse. However, it is not clear whether one can eliminate the anthropic 
argument entirely from this scenario. 

To conclnde, let us mention some important directions for fnture work. 
In this paper we have focussed on particular wrapping configurations, which 
lead to cosmologically interesting evolution. The simple topological argnment 
introdnced by Brandenberger and Vafa snffices to show that onr wrapping 
conflgnrations are possible: they conld arise from thermal flnctnations in the 
early nniverse. However we have not attempted to analyze the likelihood of 
such a fluctuation. We are currently stndying the thermodynamics of the 
early nniverse in this scenario, including the thermal creation and annihila¬ 
tion of branes |^. An important ingredient in this is calcnlating the growth 


of the cross-section with wrapped radii and temperatnre; this is analogons 
to the growth of transverse fluctuations of a string with oscillation number 
and resolntion time |^. This will allow us to estimate the probability of 


obtaining the brane wrapping conflgnrations we have investigated. 

In the present paper we have only considered spatially homogeneous uni¬ 
verses. The assnmption of spatial homogeneity greatly simplifies the analysis, 
bnt to trnly nnderstand the robustness of this scenario, it is important to 
relax this assnmption. Several new effects must be included in spatially in¬ 
homogeneous universes. In particnlar the long-range forces between between 
branes, mediated by the exchange of snpergravity particles, must be taken 
into acconnt. 
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